Boundary and Bulk Phase Transitions in the Two Dimensional Q > 4 State Potts 

Model 
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The surface and bulk properties of the two-dimensional Q > 4 state Potts model in the vicinity of 
the first order bulk transition point have been studied by exact calculations and by density matrix 
renormalization group techniques. For the surface transition the complete analytical solution of 
the problem is presented in the Q — > oo limit, including the critical and tricritical exponents, 
magnetization profiles and scaling functions. We have shown in an accurate numerical study that 
the universality class of the surface transition is independent of the value of Q > 4. For the bulk 
transition we have numerically calculated the latent heat and the magnetization discontinuity and 
we have shown that the correlation lengths in the ordered and disordered phases are identical at the 
transition point. 
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I. INTRODUCTION 

The two-dimensional Q-state Potts modelS is an im- 
portant testing ground of statistical physics for different 
approximate theories and numerical methods. For Q < 4 
the model has a second order phase transition and the 
critical exponents, both in the bulk and at the surface, 
are Q-dependent as they have been predicted by confor- 
mal invarianceo and by Coulomb-gas methodsB. 

On the other hand in the-first order regime (i.e. Q > 4) 
there are analytical results] about the properties of the 
system at the transition point, such as the latent heat, 
AE, the discontinuity of the magnetizationd, r&m, and 
the correlation length in the disordered phaseO, but 
still some issues are unsolved. One of these is the re- 
lation between the correlation lengths in the order (£ D ) 
and in the disordered (£<j) phase for which different the- 
oretical approaches which predict eitherQ £ /£d = 1/2 or 
£,o/£,d — 1, where the latter is supported by recent Monte 
Carlo simulations^!. 

Another unclarified problem of the Q > 4 model con- 
cerns the behavior of the surfaca,phase transition. In 
studies of surface critical behaviom most of the theoreti- 
cal investigations have been focused on the so-called ordi- 
nary transition which occurs at a bulk continuous tran- 
sition in absence of enhanced surface couplings and/or 
surface fields. In this case scaling ideas and, in two di- 
mensions, conformal invarianceo are often used to predict 
surface critical exponents and scaling functions. 

A completely different approach is necessary for sys- 
tems with a discontinuous bulk phase transition, like the 
Q > 4 state Potts model, where scaling and, conformal 
invariance do not hold. It was LipowskyM who first 
showed, in a mean-field calculation, that one can have 



a first order transition in the bulk and a second order 
one at the surface. This phenomenon, known as surface 
induced disorder_J.SH)) and closely related to the wet- 
ting i^henomenaliil, has been studied at the mean field 
leveMa, by means of Mont^,Carlo simulationsE-3 and in 
some experimental systcmstj, but no exact results are 
available on microscopic models so far. 

In this paper we present exact and numerical results on 
surface phase transitions in the two-dimensional Q > 4 
state Potts model. Our study is based on two comple- 
mentary approaches. First, we show that there exists 
a systematic large Q expansion of the problem, which 
is exactly soluble in the Q — > oo limit. We show that in 
this limit there is SID with a surface magnetization which 
vanishes continuously by approaching the bulk transition 
temperature and we find the exact values of the whole 
sets of surface exponents. At a sufficiently large value of 
a surface field applied to the spins at the free boundary, 
the surface transition becomes first order. The continu- 
ous and first order regimes are separated by a tricritical 
transition whose exponents are also found exactly. 

The analytical results are complemented by numerical 
investigations based on depaity matrix renormalization 
group (DMRG) techniquesll3, through which we calcu- 
late magnetization and energy density profiles in L x oo 
strips for 5 < Q < 9. Following the approach of Ref. 
|l"7| we restricted ourselves to moderately large values of 
L; typically we considered L < 40 and occasionally also 
larger systems (up to L = 80). We analyzed the numeri- 
cally accurate data by efficient extrapolation techniques 
and using information of scaling functions obtained in 
the exact solution at Q — * oo. We found again SID and 
strong indications that the surface exponents do not de- 
pend on the value of Q. 

The same type of DMRG technique is used also to cal- 
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culate properties of the bulk transition. Here we note 
that for weak first-order transitions, when Q is, close to 
four, it is a challenge for numerical methodstS (Monte 
Carlo simulation, renormalization group (RG) methods, 
finite-size scaling, series expansion, etc.) to detect the 
correct order of the transition and to estimate the corre- 
sponding discontinuities. In all cases studied we found 
good agreement between the numerically extrapolated 
values and the exact results, which once again proves 
the powerfulness and reliability of the DMRG algorithm. 
The numerical calculation of the bulk correlation length 
in the ordered phase (£d) for Q = 9 confirms with good 
accuracy the conjectured relation ^ = £ G . 

The paper is organized as follows. In Section || we 
introduce the model, define the basic quantities and out- 
line the use of the DMRG method for the problem. In 
Section III we present some numerical results concerning 
bulk quantities and analyze the behavior of the bulk cor- 
relation length. The surface transition of the model is 



studied in Section IV by analytical and numerical meth- 
ods. Section ^ concludes the paper with some discussion. 



II. FORMALISM 

We consider the Q-state Potts model in a L x oo square 
lattice defined by the Hamiltonian 
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in terms of a Q component spin variable s;^ = 
1,2, ... ,Q. We use for the inverse temperature (3=1 
and consider the ferromagnetic model K\, > 0, where 
the transition temperature is given from a duality rela- 
tion as: 
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We work in the transfer matrix formalism, with the trans- 
fer matrix along the strip given bytij: 



x exp 



T 



L—l 



exp 



L Q-l 



E E M ? 

1=1 q=l 



K 2 ( 8(8, - e l+1 ) - i 



i=i 



(3) 



where M; is a spin-flip operator: M/|s ; ) = |s ; + l,mod Q) 
and K\ is the dual coupling defined as (e Kl — l)(e Kl — 
1) = Q, thus at the critical point K* — K2- 

Different physical quantities of the model can be ex- 
pressed by the leading and next to leading eigenvalues 
of T, Ao and Ai, respectively, and by the corresponding 



dominant eigenvector \vq). For example one obtains for 
the magnetization profile 



mi 



Q (v \S( Sl - l)\v ) - 1 
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Q-l 

and for the energy density 

= (vo\$(sL k - s l+ltk )\v ) . (5) 

The correlation length along the strip is given by 

e,| =log- 1 (Ao/A 1 ) . (6) 

In what follows we use symmetry breaking boundary con- 
ditions (BC) and fix spins along one of the surfaces to 
the same state s±.k = 1, k = 0, ±1, . . . and at the other 
surface of the strip we apply a surface field of the form 

-HsE k [8(sL,k-l)-l/Q}. 

In the transfer matrix approach one needs to diagonal- 
izc the transfer matrix, which is technically very demand- 
ing, since the dimension of the matrix being Q L_1 x Q L_1 
grows very rapidly with the width of the strip. With the 
traditional finite size scaling methods based on the use 
of the Lanczos algorithm one can only treat very limited 
range of the values of L, especially for not too small Q- 
s. This limitation is the reason why density profiles and 
surface quantities have not been studied following this 
method. Another possibility is the representation of the 
Potts variables in terms of Temperly-Lieb (TL) opera- 
tors. This has the advantage that the transfer matrix 
dimension does not depend on Q, which enters in the 
calculation as a parameter, therefore this representation 
could be suitable for studying large Q. However quanti- 
ties as density profiles cannot be directly expressed by 
the TL operators. For the study of these quantities one 
should then use the DMRG in the original Potts repre- 
sentation, j— | 

In a previous papertll, hereafter referred to as Paper 
I, we have used the DMRG method for the Q < 4 Potts 
models to study the critical properties, such as density 
profiles, Casimir amplitudes and critical exponents. Here 
we extend the investigations to the first order transition 
regime with Q > 4. 

The principle and the technical details of the. DMRG 
method are well documented in the literatureo and its 
application to the Q-state Potts model is described in I. 
Here we just mention that in the DMRG method a strip 
of large width L is described by an approximated transfer 
matrix as represented in Fig. 0(a), where the rectangles 
denotes blocks containing a large number of spins, but 
described by only m states. For the Potts model the 
total dimension of the matrix would be m?Q 2 x m 2 Q 2 . 
Although m, the number of states kept, can be small 
(typically m ~ 100) the algorithm provides accurate re- 
sults even for large L. For a strip of given width, as 
Q increases, one needs to increase m in order to have 
sufficient accuracy. This limits the range of values of Q 
that can be investigated with this method. Therefore 
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for the Q = 9 model we have implemented a method, 
which in principle can be used for Q = 2 2 , 3 2 , 4 2 , . . ., 
and which considerably speeds up the numerical compu- 
tation. It consists in rewriting the original Hamiltonian 
for the Q = g 2 -Potts model as a q-Potts model with pla- 
quette and further neighbor interactions as illustrated in 
Fig. 0(b). Two neighboring g-spins (or pseudospins) in a 
row form a Q-spin; in terms of these new variables which 
we indicate as Sjj the Hamiltonian (Q) takes the following 
form: 

oo L — 1 

- f3H = Ki ^2 $(hi-i,k - s2i+2,k)S(s 2 Lk - hi+i,k) 

k— — oo l—l 
oo L 

+ K 2 ^2 y^^jsi.k - h,k+i)S(hi,k - hi,k+i) ■ (7) 

k— — oo £=1 

When the DMRG algorithm is applied to this Hamil- 
tonian one deals with matrices of dimension m 2 Q x m 2 Q, 
and the number of states kept in the blocks can be en- 
larged, with respects of the original representation. Com- 
parisons of the standard and pseudo-spin version for the 
Q = 9 case indicate that the latter provides a remarkable 
increase in computational speed and accuracy. We also 
point out that a similar idea has been applied recently to 
study otquantum systems with a large number of states 
per siteci. 
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up to a constant term and h = KI/K2Q. Thus at the 
critical point h c = 1/Q. 

In this Hamiltonian version of the model quantities of 
physical interest are derived from the ground state wave- 
function l^o) and from the ground state and first ex- 
cited state energies En and E%. In the anisotropic limit 
the transfer matrix and the Hamilton operator commute, 
[T, TL] = 0; the corresponding eigenvectors are the same, 
|*o) = I wo), consequently the density profiles in Eqs.ffl) 
and (§) stay in the same form. On the other hand the 
correlation length along the strip is now expressed by the 
inverse gap 
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(9) 



instead of Eq.(||). In the Hamiltonian version the sym- 
metry breaking BC are equivalent to fix one of the sur- 
face spins to the state \s\) — 11) and to add a term 
-H s [S(s L - 1) - 1/Q] to H in (|) in order to study the 
effect of a surface field H s on the system. This Hamilto- 
nian version of the model will be used to solve the surface 
properties of the model in the Q — ► 00 limit. 




(b) 






FIG. 1. (a) Schematic view of a transfer matrix element as 
generated by the DMRG algorithm; the blocks describe ap- 
proximately a large number of spins using m states only, (b) 
Pseudospin version of the Q — q 2 Potts model as described in 
the text; two q spins grouped by the dashed ellipses form a Q 
spin and interactions are denned on plaquettes and between 
further neighbors along the row as indicated in the figure. 

Before closing this section we turn back to the transfer 
matrix in (g|), which in the vicinity of the critical point 
can be simplified in the anisotropic or Hamiltonian limit 
of strong vertical and weak horizontal couplings where 
both K\ and K 2 are small. Then one can combine the 
exponentials to obtain T = cxp(— K2H) with a Hamilton 



III. BULK TRANSITION 

In this Section we consider the isotropic model, K± = 
K2, and present numerical results about the bulk tran- 
sition obtained by the DMRG method for 5 < Q < 9. 
Our aim in calculating the latent heat and the magne- 
tization discontinuity is to demonstrate the accuracy of 
the DMRG approach. On the other hand we used the 
numerical results about the correlation length in the or- 
dered phase at the transition point to analyze the ratio 

We start with the calculation of the latent heat, AE, 
which is just the difference in the energy densities in the 
ordered and the disordered phases at the transition point 
AE = e° — e d . The arithmetic average of the above two 
energy densities is the energy per bond, e = (e° + £ d )/f 
which is known from self-duality at the transition point! 



1 I 1 



(10) 



Consequently it is enough to calculate the energy density 
in the ordered phase and the latent heat then follows as: 



AE = 2(e° - e) 



(11) 



To estimate e° we have considered L x 00 systems with 
fixed spin BC at both edge of the strip and calculated 
by the DMRG method the energy density in (0) at the 
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middle of the strip, e/,/2, for a series of the lengths 
L = 8, 12, 16, . . . , L max . We selected moderately large 
finite systems, L max ~ 40 — 80, so that the DMRG data 
are still accurate enough (with a relative accuracy of 
10~ 5 — 10~ 6 ) to use sophisticated series extrapolation pro- 
cedures, such as the BST algorithmEl As we have shown 
in paper I, this type of combination of the DMRG method 
and the finite-size extrapolation methods results in accu- 
rate data in the thermodynamic limit. (In a very recent 
corner transfer matrix DMRG calculation of the latent 
heat of the Q = 5 Potts model Nishino and OkunishiE£l 
used another strategy and performed the calculation for 
several thousand of layers without using series extrapo- 
lation algorithms.) To demonstrate the accuracy of our 
approach we show in Table Q the BST extrapolants for 
the latent heat of the Q = 8 model. As seen in this 
Table the extrapolants seem to converge to a value of 
AE(Q = 8) = 0.242(2), which is fairly close to the exact 
value 0.2432. Results of the extrapolated latent heats 
for Q = 5, 6, 8 and 9, together with the exact values are 
presented in Table O. It is seen that the DMRG method 
combined with series extrapolation techniques predicts 
the right order of the transition and the estimated latent 
heats are in good agreement with the exact results, even 
for Q = 5. As expected, the accuracy of the numeri- 
cal results is lower for weakly first order transitions, i.e. 
which are characterized by a small latent heat. 

Next we turn to calculate the magnetization disconti- 
nuity, Am, at the transition point. Since the magnetiza- 
tion in the disordered phase is zero Am is just the bulk 
magnetization in the ordered phase. This can be esti- 
mated from the magnetization profile in (^) in the mid- 
dle of the strip, m^ 2 , with fixed spin BC at both edges 
of the strip. As for the latent heat we have calculated 
mi/2 for a series of moderately large L-s by the DMRG 
method and then performed an extrapolation with the 
BST procedure. The estimated values of Am and the 
corresponding exact values are collected in Table |j| for 
different values of Q. As seen in the Table the DMRG 
results agree well with the exact results and the relative 
accuracy in this case is even higher than for the latent 
heat. This is due to the fact that the magnetization dis- 
continuity is relatively large, even for Q = 5 Am is about 
50% of the saturation value. Considering the finite size 
corrections it is interesting to note that while these cor- 
rections for Q = 5 are approximately logarithmic, for 
Q = 9 we have a L~ 2 correction behavior. 

Finally, at the bulk transition point we have calculated 
the correlation length in the ordered phase, £ Q , for Q = 9, 
since for smaller Q-s the correlation length exceeds the 
available size of the strip. 

Figure || shows a plot of the correlation length t; (L) 
for Q = 2, 3, 5 and 9, calculated from the ratio of the 
largest and the second largest eigenvalues of the transfer 
matrix in Eq. (^J). The Q = 2 and 3 cases follow the. 
expected behavior derived from conformal invariancetl 
for a strip with fixed and equal spins at the boundary, 
i.e. £ (L) — L/(2ir). Notice that also in the Q = 5 



case £ (L) follows a straight line with the same slope as 
for the Q = 2 and 3. The transition, as pointed out 
above, is only weakly first order in the case Q = 5 with 
a correlation length much larger than the largest sys- 
tem size available here. The data for Q — 9 show a 
clear sign of saturation for large L. The extrapolation 
to L — > oo has been performed by the BST method and 
yields £ Q = 15.0(4). This quantity is in good agreement 
with the exactly known correlation length of the disor- 
dered phase which is £d = 14.90. These results corrobo- 
rates the recent conjecture of Ref. |^, namely that £ = £,d, 
based on results from Monte Carlo simulations. 
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FIG. 2. Correlation length £a(L) for fixed/fixed boundary 
conditions as function of the strip width L. The dotted line 
is the conformal result. 



IV. SURFACE TRANSITION 

The surface transition of the model is studied by two 
complementary methods. First, we present the analyti- 
cal solution of the problem in the Q — ► oo limit, then the 
transition for finite Q-s are studied numerically by the 
DMRG method. Also in this numerical part of the work 
we use the analytical results obtained for Q — > oo about 
the critical singularities and scaling functions, as a help 
in analyzing the DMRG data. 



A. Analytical solution in the Q — > oo limit 

Here we work with the anisotropic version of the model, 
which is characterised by the Hamiltonian in (|^) and use 
the BC described below (||), i.e. fixed spin at i = 1 and a 
surface field of strength H s in the other end, i = L. (For 
a free surface we have H s = 0.) 

The Hamiltonian TL is represented by a matrix of di- 
mensions x which can be further reduced in 
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a block-diagonal form by making use of the symmetries 
of the problem. The dimension of the block which con- 
tains the ground state (and also -the first excited state) 
is independent of the value of CjE3. For small systems, 
L = 2 and L = 3, this block, which has the highest sym- 
metry, has been explicitly constructed in the Appendix. 
Analyzing the structure of this ground-state block one 
can notice further simplifications in the vicinity of the 
critical point, where one can perform a large Q expan- 
sion, such that the corrections to the leading behaviour 
are of 0(Q^ 1 ^ 2 ). As shown in the Appendix in the limit 
Q — ► oo the ground-state block has a dimension of L with 
an eigenvalue matrix H given by 



H = — (L - 1) - hy/Q 



\0 



1 

2t 



\ 



1 (i-l)t/ 



(12) 

Here t = (hQ — l)/hy/Q measures the distance from the 
bulk transition point and h s = H s /hy/Q. Denoting the 
components of the eigenvectors of H as ipi (n) with n = 1, 
. . . , L one can express the magnetization profile in terms 
of the ground state wavefunction as follows: 



L + l-l 



mi 



(13) 



n=l 



At the transition point t = the matrix (12) is exactly 
diagonalizable with eigenvectors 

ipi(ri) = AiSm[ki(L + 1 - n)}, (14) 

and eigenvalues 

Ei = -(L - 1) - 2hy/Qcosk l , (15) 

where i = 0, 1, . . . , L — 1 and the wavenumbers fc$ are 
determined by the equation: 



h s sm(Lki) = sin[(i + l)ki 



(16) 



As known by standard analysis the spectrum and the 
ground state of H have different properties for h s < 1 
and h s > 1, respectively. In the following we treat the 
different regions of surface transition separately. 



1. Critical surface transition regime - h s 



< 1 



L + l 



(17) 



and the normalization in ((lj) is A l = [2/(L + 1)] 1/2 . 
The boundary magnetization from (jl 3|) is given by: 



m s = m L 



. 2 

■ sin 



L+l \L+1 



(18) 



thus it goes zero at the bulk transition point, t = 0. Con- 
sequently the surface phase transition is of second order, 
furthermore the anomalous dimension of the surface mag- 
netization operator is 



(19) 



as obtained from the finite size dependence in (|l^ 

The magnetization profile from (^) can also be eval- 
uated, which in the continuum approximation L,! > 1 
and as a function of the distance from the free surface, 
Al = L + 1 — I, is given by 



Al 1 f Al 
rriAi = — — sm | 27r— 

L Z7T 



AIY 



This profile is shown as a solid line in Fig. |^. 
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For h s < 1 there are L real solutions of Eq. (|T^) and the 
spectrum of H is quasi-continuous. The qualitative be- 
haviour of the surface transition is the same in the whole 
region, therefore we consider here the free surface case, 
h s = 0, when the the solution of ( |l6| ) is given explicitly 
as: 



FIG. 3. Magnetization profiles of the Q — * oo Potts model 
at the bulk transition point with fixed spin BC at I = 1 and 
with different values of the surface field, h s , at I = L. Solid 
line - critical surface transition, h s — 0; dotted line - tricriti- 
cal surface transition, h 3 — 1; dashed line - first order surface 
transition, h s > 1. 
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In a strip of width L the correlation length parallel 
with the free surface is given from (^) : 



C7 1 = 2hy/Q 



cos 



L + l 



2tt 



(21) 



Since the correlation length perpendicular to the strip is 
just £j_ ~ L, we have the asymptotic relation £|| ~ , 
with an anisotropy exponent 



= 2 



(22) 



Away from the bulk transition point for < —t ^ 1 
the boundary magnetization can be studied by a pertur- 
bation expansion, when the diagonal term of H in (^), 
, j = —t(i—l)6ij, is considered as a perturbation. From 
this calculation the first order correction to the boundary 
magnetization is given for L — > oo and < — t <C 1 as: 



2^ 2 1 /7T 



(23) 



Thus the surface magnetization exponent, defined in the 



thermodynamic limit [L 
by: 



oo) as m s ~ {—t) s is given 



(24) 



The small i-behavior of the system can also be investi- 
gated in the frame of a continuum approximation, when 
the eigenvalue problem of H in (|l2"|), which is a second- 
order difference equation, is transformed into a differen- 
tial equation of the continuous variable x = i: 



dx 2 



+ (e + tx)tp(x) = , 



(25) 



with e = 2 + E/{h\fQ) and t — t/(h^/Q), whereas the 
boundary conditions are ip(x = 0) = iMx = L) = 0. 
From the properties of the solution of (g5|) one obtains 
for the temperature dependence of the correlation lengths 
as: 



Cii ~ \t\ 



-2/3 



-1/3 



(26) 



thus the corresponding critical exponents are i^" = 2/3 
and i/j_ r = 1/3 and the scaling relation v\\/v±_ = z CI is 
satisfied with z cr in (p2|). The exponents of the surface 
critical transition are collected in Table III. 



2. First-order surface transition regime - h s > 1 

The qualitative behavior of the solution of H is chang- 
ing when the strength of the surface field, h s grows over 
h s = 1. For /i s > 1 the ground state of H has the form 
of a localised mode and the corresponding wavenumbcr 



in (16) is complex, fco = in. In the limit L — > oo the 
wavefunction of this localised mode is given by 



^ loc (n) = (h 2 s lf/ 2 h- n 



(27) 



with the energy 

E loa = -(L-l)-hy/Q(h s +h- 1 ) , (28) 

which is separated from the continuum part of the spec- 
trum given in the form of (15). 

The surface magnetization at the bulk transition point, 
t = 0, from (||) 



ms = 1 -h 2 . 



(29) 



is finite, thus the surface transition is of first-order for 
h s > 1. The S^i exponent defined as m s ~ (Shs) 1 ^ 1 - 1 
with Sh s = h s — 1 is given by: 



(30) 



where the superscript "tr" refers to the tricritical point, 
which is approached as h s — > 1 + . 

The surface correlation length parallel with the surface 
is finite in the L —> oo limit 



Ci 1 =E X - E loc = h^/Q{h s + h; 1 - 2) ~ (Sh s ) 2 , (31) 



as expected for a discontinuous surface transition. In 
( |3l| ) for Ei the edge of the continuum spectrum in Jl5| ) 
is used. The magnetization profile as a function of the 
distance from the free surface, Al = L + 1 — I, is given 
by: 



rriAi = 1 - h 



-2AZ 



h. > 1, L 



(32) 



As seen from the dashed line of Fig. || there is a surface 
region of the profile having a size of 



i 



21n/i, 



(SKy 1 



(33) 



which defines the perpendicular correlation length at the 
surface. As the tricritical point is approached, h s — * 1 + , 
the two correlation lengths are related as £^ ~ ^m, thus 
the anisotropy exponent in the tricritical point 



z tr = 2 



(34) 



is the same as in the critical region Eq.(G2T). 
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3. Tricritical surface transition - h s = 1 



4- Crossover regime 



As we mentioned before the tricritical point with h s = 
1 separates the continuous and first-order surface tran- 
sition regimes. At the bulk transition point, t = 0, the 
solution of (Hq) for the wavenumbers is given by 



2L+ 1 



(2* + l) 



1 , 



(35) 



and the normalization in ( |l4| ) is Aj = 2/ v / 2L + 1. 
The surface magnetization 



2L+ 1 



2(2Z + 1) 



(36) 



goes to zero as L — ► oo, however the tricritical surface 
magnetization scaling dimension 



x'l = 1 , 



(37) 



is different from its value in the critical region (|l9|). 

The tricritical magnetization profile can be evaluated 
from (0), which in the continuum approximation is given 
by 



m A i = 



Al 


i 


T " 


h - 

7T 


Al 








T " 



A? 



(38) 



as a function of the distance from the free surface, Al. 
This profile is also shown in Fig. |as a dotted line where 
one can notice the different behavior close to the free 
surface for the critical and tricritical cases. This is due 
to the fact that the profiles grow with the respective x s 
scaling dimensions, which are different in the two fixed 
points. The correlation length parallel with the surface 
can be calculated as in (|2l| ) with the asymptotic behav- 
ior £|| = L 2 . Thus the anisotropy exponent is z tr = 2 in 
accordance with the previous derivation in (j34|). 

Away from the bulk transition point for \t\ <C 1 one 
can repeat the calculations as described for the critical 
surface transition. In the continuum approximation the 
differential equation in ( p5| ) remains the same and just 
the boundary condition at x — will be changed to 
dip/dx\ x= o = 0. The tricritical exponents obtained from 
the solution of the differential equation are presented in 
Table III, together with their critical counterparts. We 
stress that scaling in the surface region is anisotropic, 
both in the critical and tricritical fixed points the cor- 
relation lengths are related as £|| ~ with the same 
anisotropy exponent z = 2. On the other hand the x s and 
(3 S exponents are different in the two cases, but the scal- 
ing relation S — u±x s is satisfied in both fixed points. 
Notice also that the exponents in Table III are the same 
as for the restricted solid-on-solid mpdcl of the interface 
localization-delocalization transitioned. 



It is interesting to study the behavior of the surface 
magnetization in the vicinity of the tricritical point when 
< 1 — h s <C 0. This is a crossover region where the sur- 
face critical behavior for finite L is strongly influenced 
by the presence of the tricritical transition occurring at 
h s = 1. 



Analyzing the exact relations in (fDl-16) one can define 
a length scale £ s = l/2|ln/i s | which is identical to £j_ 
defined above in the first order case h s > 1 in ([33|). Then 
one can show that the finite size surface magnetization at 
the transition point, t — 0, satisfies the scaling relation: 



m s (L,hih s ) = L Xs m s (L/^ s ) 



(39) 



with £ s defined above and where the scaling function 
fh s {xf) goes to a constant for y — > and behaves as 
m s (|/) ~ 1/y 2 for y — > oo. This scaling function is ex- 
pected to be universal, provided the surface critical be- 
havior of the model is independent of Q > 4. We shall 
study this point numerically in the next Section. 



B. DMRG study 



Having the complete analytical solution for Q — > oo at 
hand, we are interested in the surface critical behavior for 
finite values of Q. The expected surface phase diagram 
of the model is schematically drawn in Fig. |]. 




FIG. 4. Schematic phase diagram for the surface critical 
behavior of the Q-state Potts model as function of the ap- 
plied surface field h s . The thick solid line departing from 
Q = 4 denotes a tricritical line, which separates a continu- 
ous from a first order surface transition. For Q < 4 (along 
the thick dashed line) the surface critical exponents depend 
continuously on Q. 
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In the second-order bulk transition regime, Q < 4, the 
critical and tricritical surface transition lines merge and 
are characterized by Q-dependent exponents. For Q > 4, 
in the first-order bulk transition regime, the second- and 
first-order surface transition regions are separated by a 
tricritical line, which is expected to depart from Q = 4. 
It should also be noticed that some of the tricritical expo- 
nents obtained from the exact solution at Q — > oo are the 
same as for the Q = 4 Potts model: x* r — x s (Q = 4) and 
v \\ = t / buik(<3 = 4), whereas other are different. Then 
one can expect two scenarios about the surface (critical 
and tricritical) exponents of the Q > 4 model: i) the 
exponents vary with Q, like for Q < 4; or ii) the expo- 
nents are universal and their fixed-point values are those 
calculated at Q — > oo. 

To decide between the two possibilities we have in- 
vestigated the problem numerically using the density 
matrix renormalization group. As for the bulk tran- 
sition in Section III we have considered the isotropic 
model for 5 < Q < 9 on moderately large finite strips 
(L < 40) where we have calculated the magnetization 
profiles with fixed-free boundary conditions, using the 
finite system algorithm described inE3. Comparing the 
numerical data foij-different number of states kept in the 
DMRG procedureEj we could estimate the accuracy for 
the finite surface magnetization as about 5 — 6 digits. 

Figures |^ and ^ show plots of magnetization profiles 
for the Q = 7 and Q = 9 cases, for different values 
of the strip width L. From the surface magnetization 
m s (L) at the bulk transition point and at h s = we 
first formed finite-size extrapolants: x s (L) — ln(m s (L — 
4)/m s (L))/ln(l - A/L) from the data for L = 8, 12, 16, 
. . . 40. For L — ► oo one has x s (L) — > x s . Table IV shows 
the values of x s (L) for the cases Q = 8 and (5 = 9. All 
the finite size data are above the value x s = 1 (the sur- 
face critical dimension for the Q = 4 model) and tend 
to increase at larger Q. ExtraBolations from the BST 
sequence extrapolation methocO yield x s = 2.7(4) for 
Q = 8 and x s — 2.9(3) for Q = 9, both compatible with 
the value of x s — 3 found in the Q — > oo limit. These 
estimates have quite large error bars (from this reason 
we could not get a sensitive estimate for Q < 7) and we 
did not find nice convergence of the BST extrapolants 
for x s (L) as was found in the cases Q = 2 and Q = 3 
in Ref. [IT]. This is typically due to strong finite size cor- 
rections, which in the present case are essentially caused 
by two effects. First of all the bulk correlation length 
£b(Q), which remains finite at a first-order transition is 
still very large for the values of Q analyzed. As seen in 
the values of the bulk correlation lengths in Table || only 
the highest values of Q are reliable for extrapolations of 
finite size data with L < 40. A second point is that in 
the range 5 < Q < 9 we expect that the tricritical sur- 
face transition is at small surface fields, so that even at 
h s = strong influence from the tricritical fixed-point 
can be observed. 




FIG. 5. Magnetization profiles at the bulk transition point 
for fixed/free boundary conditions and for L — 8, 12, 16, . . . , 
40 with Q = 7. 
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FIG. 6. As in Fig. | but for Q - 



To clarify this point we have reanalyzed our numerical 
data using the analytical scaling relation in (|3^) for the 
Q — > oo limit, which describes the behavior of the sur- 
face magnetization at and in the vicinity of the tricritical 
point. As already mentioned the scaling function in (|39| ) 
is expected to be universal, provided the surface critical 
behavior of the model is independent of Q > 4. In the 
analysis of numerical results we have compared the ratio 
of surface magnetizations for a finite Q and h s = 0, with 
the corresponding analytical results at Q — > oo and at 
finite h s . We found the value of h s such that the equal- 
ity: m?(L)/m?(L-4) = mf (L, h s )/mf{L - 4, h s ) is 
satisfied. The universality hypothesis about the scaling 



8 



function is connected with the convergence of the finite 
size estimates h s (L) to a finite limit. Moreover, the dif- 
ference, Ah s = 1 — h s , gives an estimate on the position 
of the tricritical line at a given Q > 4. 

As seen in Table [V the Ah s (L) finite size extrapolants 
for the Q = 7 model are nicely converging to a limiting 
value and the same observation has been found for the 
other Q-s, as well. The extrapolated position of the tri- 
critical line is found as Ah s = 0.0125, 0.0214 0.0310 and 
0.040, for Q = 6, 7, 8 and 9, respectively^. This is a 
strong numerical indication that the surface critical phe- 
nomena of the Q > 4 state Potts model is universal and 
the corresponding exponents are those given in Tabic III . 
The influence of the tricritical point is also visible in the 
profiles for Q = 7 of Fig. || where the magnetization 
has an approximately linear behavior as function of the 
distance from the free surface as found for the tricritical 
profile in the limit Q — > oo (see dotted line of Fig. ||). 
In the case Q = 9 instead, the curvature of the profile is 
clearly noticeable for the largest systems analyzed. 

Finally, we turn to compare the two finite length scales 
appearing in the problem: the £ s ~ \/2Ah s surface 
length scale and the £b bulk correlation length. Assuming 
anisotropic scaling as ~ £j , we obtain 6 = 1.15, 1.24 
and 1.4 for Q — 8,7 and 6, respectively. In the Q — > 4 + 
limit, when the correlation lengths are divergent, it is 
natural to assume that the anisotropy exponent will ap- 
proach 8 = 2, so that the perpendicular and parallel cor- 
relation lengths are related as in the critical and tricritical 
surface transitions (see Table III). With this assumption 
we can conjecture the position of the tricritical line as 



lim AhJQ) 



MQ)] 



-1/2 



exp 



2 v / Q-4 



(40) 



where the limiting form of £b(Q) can be found ir 



V. DISCUSSION 

In this paper the two-dimensional Q-state Potts model 
has been studied by analytical and numerical methods in 
the first-order bulk transition regime, Q > 4. In the bulk 
transition point we have used the DMRG method and se- 
quence extrapolation techniques to accurately estimate 
the latent heat and the magnetization discontinuity at 
the transition point. Calculating the correlation lengths 
at the transition point we have shown with good numer- 
ical accuracy that they are identical in the ordered and 
in the disordered phases. 

For the surface transition we have demonstrated the 
existence of surface induced disorder, for the first time 
for a microscopic model. We have calculated the com- 
plete set of critical and tricritical surface exponents and 
showed that the surface transition of the model is uni- 
versal for Q > 4, thus do not depend on the value of the 
discontinuities at the bulk transition point. Based on this 
observation we conjecture that the surface exponents in 
Table III are universal for all two dimensional systems 
with a first-order bulk transition. This statement is in 



accordance with the prediction of the restricted solid- 
on-solid model of the interface localization-delocalization 
transitioned. It is interesting to note that the transfer 
matrix of the restricted solid-on-solid model is formally 
equivalent to the Hamiltonian H in ([l2|), which repre- 
sents the ground-state block of the anisotropic version of 
the Q — > oo Potts model. The physical interpretation 
of the same operator in the two problems is, however, 
completely different. 
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APPENDIX: SOLUTION OF THE 
GROUND-STATE SECTOR OF H 

Here we show how to construct the ground-state sector 
of (||) for the smallest systems, L — 2 and L = 3 and then 
use these results to define a systematic 1/Q-expansion in 
the vicinity of the bulk transition point. 

The ground-state sector for the L = 2 system is 
spanned by two vectors: 



i>2 



n 



i 



(12 + 13 + . 



1?) 



(Al) 



and the eigenvalue matrix with zero surface field, H s = 0, 
is given by: 



H 2 



-1 



-hVQ~T -h(Q-2) 



(A2) 



For L = 3 the ground-state sector is spanned by five 
vectors: 



- 111 

1p2 = 



1 



^3 = ^(121 + 131 
-04 = ,^ ^ (122 +133 



= (112 + 113 + . . . + llg) 
. + lql) 

■ + Iqq) 
(123+ 124+... + 12g 



V(Q-i)(Q-2) 
+ 132... + lg(g-l)) (A3) 

and the eigenvalue matrix for H s = is given by: 
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H 3 = 



/ 


-2 


- hy/Q - 1 


-hy/Q 


- 1 







-hy/Q - 1 


-l-ft(Q-2) 







-ft 




-hy/Q - 1 





-ft(Q- 


-2) 


-ft 







-ft 


-h 




-1 


V 





-hy/Q - 2 


-hy/Q 


- 2 


-h2y/Q 



-hy /Q~~2 
-hy/Q~2 



-h2y/Q^2 

-ft2(Q - 3) y 



(4) 



It is remarkable that the dimension of the ground- 
state sector is independent of the value of Q and only 
depends on L. In the vicinity of the bulk transition 
point h c = 1/Q further simplifications take place in 
the Q — > oo limit. To see this, first, we consider the 
L = 2 problem, where the two vectors ipi and ip2 are 
now degenerate with the same diagonal matrix-elements, 
whereas the off-diagonal matrix-elements are —hy/Q in 
leading order of Similarly, in the L = 3 sys- 

tem there are thre e a symptotically degenerate states 
ipi, "02 and "05 in (A3) with diagonal matrix-elements 
Hi! = H22 — H§ 5 = — 2 and these states are separated 



l/j 4 



11. ..1 

—^(11. ..12 + 11. ..13 



from the others by an amount of 0(1). Since the off- 
diagonal matrix-elements between the relevant states are 



H 21 



H 25 



-hy/Q 1/y/Q the two 



non-relevant states do not influence the behavior of the 
ground state (and the first excited state) in the large Q 
limit. 

Generalizing the observation found in the L = 2 and 
L = 3 systems one can show that in the vicinity of the 
bulk transition point in the Q — > 00 limit for general L 
the ground state sector is spanned by L asymptotically 
degenerate vectors: 



Q lf 

— (11. ..123 + 11. ..24 

-^-(ll... 1234+ 11... 1235 
Q i/Z 



. + 11...1Q) 

+ 11...Q(Q-1)) 

. + ll...lQ(Q-l)(Q-2)) 



0i 



Q(i-i)/ 



^(1234 . . . L + . . . + Q(Q - 1) . . . (Q - L)) 



and the corresponding eigenvalue matrix is given in (|l2]). 
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TABLE I. Table of the BST extrapolants for the latent heat for the Q — 8 Potts model. In the first column are the DMRG 
data for finite systems with L = 8, 12, 16, ... , 40. (The parameter of the BST algorithm is e = 2.) 



0.335380 


0.288286 


0.312288 


0.236212 


0.242762 


0.242690 


0.242082 


0.307476 


0.272438 


0.227422 


0.243086 


0.242684 


0.242082 


0.241738 


0.291096 


0.263010 


0.240970 


0.242728 


0.242084 


0.241736 


0.241674 


0.280320 


0.256912 


0.242214 


0.242068 


0.241752 


0.241670 




0.272726 


0.252760 


0.242126 


0.241644 


0.241686 






0.267128 


0.249838 


0.241868 


0.241434 








0.262864 


0.247734 


0.241658 










0.259538 


0.246198 












0.256894 















0.241738 
0.241674 



0.241674 



TABLE II. Physical quantities at the bulk transition point 
of the Q state Potts model as found from finite size BST ex- 
trapolation of DMRG data (* denotes the exact values). 



Q 



AE 



AE (*) 



Am 



Am w 



Co 



5 0.028(3) 

6 0.102(3) 

8 0.242(2) 

9 0.299(2) 



0.0265 
0.1007 
0.2432 
0.2998 



0.49(1) 
0.66(1) 
0.798(3) 
0.833(2) 



0.4921 
0.6652 
0.7998 
0.8333 



15.0(4) 



A*) 

Sd 



2512. 

158.9 
23.9 
14.9 



TABLE III. Surface critical and tricritical exponents of the 
Q state Potts model calculated in the Q — > oo limit. 



TABLE IV. Finite size extrapolants for the surface mag- 
netization scaling dimension, x®(L), and for the position of 
the tricritical line, Ah^(L), as a function of the width of the 
strip for various values of Q. The last row shows the BST 
extrapolations at L — > oo. 



L 




*?= 9 (£) 


Ahf= 7 (L) 


20 


1.15892328 


1.25272509 


0.013300 


24 


1.22149616 


1.32785249 


0.017501 


28 


1.27464891 


1.39259493 


0.019548 


32 


1.32114750 


1.44961396 


0.020599 


36 


1.36226612 


1.50098953 


0.021121 


40 


1.39920082 


1.54783012 


0.021344 


oo 


2.7(4) 


2.9(3) 


0.0215(1) 



critical 



tricritical 



Ps 

"II 

Z 



3 

1 
1/3 
2/3 

2 



1 

1/3 
1/3 
2/3 
2 
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